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ABSTRACT 

Equilibrium statistical mechanics of two-dimensional flows provides an explanation and a prediction 
for the self-organization of large scale coherent structures. This theory is applied in this paper 
to the description of oceanic rings and jets, in the framework of a 1.5 layer quasi-geostrophic 
model. The theory predicts the spontaneous formation of regions where the potential vorticity is 
homogenized, with strong and localized jets at their interface. Mesoscalc rings are shown to be 
close to a statistical equilibrium: the theory accounts for their shape, their drift, and their ubiquity 
in the ocean, independently of the underlying generation mechanism. At basin scale, inertial states 
presenting mid basin eastward jets (and then different from the classical Fofonoff solution) are 
described as marginally unstable states. In that case, considering a purely inertial limit is a first 
step toward more comprehensive out of equilibrium studies that would take into account other 
essential aspects, such as wind forcing. 



1. Introduction 

Large scale coherent structures are ubiquitous in the 
ocean. Understanding the physical mechanism underlying 
their formation and persistence remains a major theoretical 
challenge. 

At mesoscalc, oceanic turbulence is mostly organized 
into westwar d propagating rings , as for instance revealed 
by altimetry (jChelton et al. l2007t) . Since typical eddy turnover 
times are much shorter than dissipation and forcing time 
scales, these rings can be studied in the inertial limit, for 
which forcing and dissipation are neglected. 

At basin scale, the dynamics are strongly influenced 
by forcing and dissipation: wind forcing plays the leading 
role in setting the gyre structures through the Sverdrup 
balance, and the concomitant effect of planetary vorticity 
gradie nts and dissipat ion explains their westward intensifi- 
cation (|Pedloskvl [l998') . Because none of these mechanisms 
are conservative processes, the inertial approach does not 
take these essential aspects into account. Conversely, ex- 
isting theories give no clear explanation of the existence 
of strong and robust eastward jets in the inertial part of 
these currents. The classical wind driven ocean theory and 
the inertial approach both give an incomplete picture, and 
complement each other. A useful step towards a compre- 
hensive non-equilibrium theory that would combine both 
approaches is to study mid-basin eastward jets in the iner- 
tial limit. Such is the focus of this paper. 



On the one hand, the problem of the self-organization 
of a turbulent flow involves a huge number of degrees of 
freedom coupled together via complex non-linear interac- 
tions. This situation makes any deterministic approach 
illusory, if not impossible. On the other hand, there can be 
abrupt and drastic changes in the large scale flow structure 
when varying a single parameter such as the energy of the 
flow. It is then appealing to study this problem with a 
statistical mechanics approach, which reduces the problem 
of large-scale organization of the flow to the study of states 
depending on a few key parameters only. 

Such a theory exists: this is the Robert-Sommcria- 
Miller (RSM hereafter) equilibrium statistical mechanics 



( Robertlll99(ll : lMillerll99nl : lRobertlll991 : lob crt and Sommerial 
199ll) . From the knowledge of the energy and the global 
distribution of potential vorticity levels provided by an ini- 
tial condition, the theory predicts the large scale flow as 
the most probable outcome of turbulent mixing. 

Here we ask the following question: can rings and jets 
be interpreted as RSM statistical equilibria in the frame- 
work of an 1.5 layer quasi-geostrophic (QG) model? 

The first attempt to use equilibrium statistical mechan- 
ics ideas to explai n the self - organ ization of 2D turbulence 
was performed by lOnsagerl (|l949h in the framework of the 
point vortex model. In order to treat flows with continu- 
ous vorticity fields, another approach has been proposed 
by Kraichnan in the framework of the truncated Euler 
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equations ( Kraichnan and Montgomery 198C)() . which has 
then been applied to quasi-geostrophic flows over tqpogra - 
phy ( Salmon et al. 19761 Carnevale and Frederiksen 1987 ). 



The truncation has a drastic consequence: only the energy 
and the enstrophy are conserved quantities, while any func- 
tion of the vorticity is conserved for the Euler equation. 
The energy-enstrophy statistical theory predicts the emer- 
gence of large scale mean flows above topography, char- 
acterized by a linear relationship between streamfunction 
and potential vorticity. 

The e xisten ce of such a linear relation was assumed by 
Fofonofll (| 1954T ) for analytical convenience, in earlier work 
on inertial ocean circulation, independently of statistical 
mechanics approaches. Fofonoff was able to compute ex- 
plicitly an inertial solution in the low energy limit. The 
emergence of such flows has then been observed in numer- 
ical simulations of freely evolving quasi-geostr ophic flows 
(|Zou and Hollowavll994tlWang and Vallisll 19941 ). The energy- 
enstrophy statistical theory has therefore been proven suc- 
cessful to interpret these Fofonoff flows as statistical equi- 
libria. However, Fofonoff flows are not observed in the real 
ocean. 

There exists actually a richer variety of energy-enstrophy 
states, as revealed by the co mputations of solu ti ons i n 
various configuration, see e.g. Mai da and Wang ( 2006h : 
Frederiksen and O'Kanel ( 2008 ) and references therein. The 
computation of any equilibrium state characterized by a 
linear q-psi relation for a given energy and circulation (which 
includes the previous solutions, but not only) has been per- 
formed recently for a l arge class of models includ i ng th e 
1.5 layer QG equations (jVenaille and Bouchedl200flll20inh . 
None of the observed inertial features of oceanic flows, such 
as coherent rings and mid-basin eastward jets, were found 
in this class of equilibrium states. 



e.g. Chavanisl ( 2008f ). Actually the equilibrium statistical 
mechanics predicts the ratio of potential vorticity invari- 
ants that remains to the large scale and the ratio that cas- 
cade to smaller and smaller scales for an inviscid dynamics. 

Computing the RSM equilibrium states requires the res- 
olution of a variational problem with an infinite number of 
constraints. This practical difficulty has been overcome by 
treating canonically other invariants than the en ergy, which 



leads to a variational problem easier to solve (jEllis et al 
200(1) . This approach is sometimes referred t o as statisti- 



cal mechanics with prior vorticity distribution (ITurkington 
19991 lEllis et all [20021: iMaida and Wand 120061: ~~ 



Chavanis 



2008). Any equilibrium state with prior vorticity distri- 



bution can be interpreted as an RSM equilibrium state 
( BouchedEioil ). This is therefore a very useful trick to 
compute the equilibrium states. However, we think that 
a physical interpretation of the prior vorticity distribution 
is problematic, since it would require to define what is a 
bath of potential vorticity. In the case of an isolated system 
(for instance a freely evolving inviscid flow), the relevant 
ensemble to consider is the one taking into account all the 
constraints of the dynamics. For this reason, we ultimately 
interpret the equilibrium states i n term of the RSM theory. 

Several studies, among wh ich lAbramov and Maidal(<2003l) : 
Dubinkina and Frank! (|2010l ). have specifically addressed 



the importance of higher potential vorticity moments for 
the equilibrium states. Taking into account these addi- 
tional dynamical invariants provides a much richer vari- 
ety of equilibrium states than previous energy-enstrophy 
theories. It has been proven useful to describe t he strato- 
spheric polar vortex (jPrieto and Schuberd 12001), the self- 

organ ization following deep convec tion events dPibattista and Majda 
20001) , or Jupiter Red Great Spot (jBouchet and Sommeria 



2002t ITurkington et al.ll200lD . 



For the Euler or QG dynamics, a major drawback of 
energy-enstrophy statistical theories is the loss of the ad- 
ditional invariants of the dynamics. A first attempt to in- 
clude the effect of higher o r der in variant was proposed by 
Carnevale and Frederiksen! (|l987l) . The generalization of 
Onsager's ideas to the Euler and quasi-geostrophic equa- 
tions with continuous vorticity field, taking into account 
all invariants, has lead to the RSM theory (IR obcr l | l990l : 
illerlll990HRoberdll99ltlRobert and Sommerial|l99ll). see 
Evink and Sreenivasanl (|2006l ); lMaida and Wand (|2006l) 
( 201lh 'or recent reviews and references. 



Marst 



An essential point of the RSM approach is that it makes 
a distinction between fine-grained potential vorticity distri- 
bution on one hand, which is conserved by inviscid flows, 
and coarse-grained distribution of potential vorticity on the 
other hand, which is not conserved. This is in complete 
agreement with the classical and well documented fact that 
potential vorticity invariants cascade towards smaller and 
smaller scales (both for inviscid and dissipative flows), see 



Most of oceanic coherent structures are surface inten- 
sified, with most of their kinetic energy located above the 
thermoclinc. In addition, eastward jets and rings are char- 
acterized by a jet width of the order of the first baro- 
clinic Rossby radius of deformation. In this paper, we con- 
sider the simplest ocean model that takes into account this 
vertical structure and this typical horizontal length scale, 
namely an equivalent barotropic, 1.5 layer QG model. 

We consider the limit of small Rossby radius of defor- 
mation, which allows analytical co mputations of statistical 
equilib ria, following the work of iBouchet and Sommeria 
( 20021 ). This assumption provides important insights for 
more general situations, even when such a scale separation 
does not exist. This is also a first step before considering 
the shallow water model, which is consistent with the scale 
separation between the Rossby radius of deformation and 
the domain scale. 

In the limit of small Rossby radius of deformation, it 
has been shown that the computation of RSM statistical 
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equilibria can be simplified i nto a Van-der -Waals Cahn- 
Hilliard variational problem (jBouchetl 120081 ). These vari- 
ational problems explain the formation and the shape of 
bubbles in thermodynamics. The existence of this formal 
analogy has been very fruitful in the description of Jovian 
vortices. This paper puts forward this approach in the 
oceanic context. 

The paper is organized as follows. Equilibrium statisti- 
cal mechanics of the 1.5 layer QG model is presented in the 
second section. The method to compute analytically statis- 
tical equilibrium states in the limit of small Rossby radius 
of deformation is presented on the third section. It allows 
for a justi fication of the potential vo rticity homogenization 
theory of iRhines and Youngl (|1982l ) without invoking any 
dissipation mechanism. The application to oceanic rings 
is discussed in a fourth section, by considering the case 
of a zonal channel on a beta plane. The application to 
mid-basin eastward jets is discussed in a fifth section, by 
considering the case of a closed domain. Notations and 
symbols are referenced in table [T] 

2. Statistical mechanics of the 1.5 layer, equivalent 
barotropic QG model 

a. The 1.5 layer QG model, its dynamical invariants and its 
dynamical equilibria 

The simplest possible inertial mid-latitude ocean model 
taking into account the stratification of the oceans and the 
sphericity of the Earth is considered in this paper. This is 
the unforced, undissipated, 1.5 layer QG model on a beta 
plane: 



dq 
dt 



v ■ Vg = 0, with v 



and q = V 2 </> 

K 



x Vip 



(1) 



(2) 



For the boundary conditions, two cases will be distinguished, 
depending on the domain geometry T>. In the case of a 
closed domain, there is an impermeability constraint (no 
flow across the boundary), which amounts to a constant 
streamfunction along the boundary. To simplify the pre- 
sentation, the condition ip = at boundaries will be consid- 
ered3- In the case of a zonal channel, the streamfunction ip 
is periodic in the x direction, and the impermeability con- 
straint applies on northern and southern boundaries. In 
the remaining, length scales are nondimcnsionalizcd such 
that the domain area \T>\ is equal to one. 



lr The physically relevant boundary condition should be ip = ipf r 
where ipj r is determined by using the mass conservation constraint 
f dr ip = {ip is proportional to interface variations) . Taking ip = 
does not change quantitatively the solutions in the domain bulk, but 
only the strength of boundary jets. 



According to Noether's Theorem, each symmetry of the 
system is associ ated with the e xistence of a dynamical in- 
variant, see e.g. Salmon! ( 1998} ). These invariants are cru- 
cial quantities, because they provide strong constraints for 
the flow evolution. Starting from ([TJ , ([2]) and the aforemen- 
tioned boundary conditions one can prove that QG flows 
conserve the energy: 



E = - 



dr 



v 



0y)tl>. (3) 



Additionally, the QG dynamics (JTJ is a transport by 
an incompressible flow, so that the area 7 (er) da occupied 
by a given vorticity level a is a dynamical invariant. The 
quantity 7(0") will be referred to as the global distribution 
of potential vorticity. The conservation of the distribution 
7 (a) is equivalent to the conservation of any moment of 
the potential vorticit y f v dr q n , a nd is rel a ted to particle 
relabelling symmetry Ripa ( 1981 ): Salmon ( 19981) . 



The stationary points of the QG equations ([T]), referred 
to as dynamical equilibrium states, satisfy v • Vg = ViJj x 
Vg = 0. It means that dynamical equilibria arc flows for 
which streamlines are isolines of potential vorticity. Then, 
any state characterized by a q — ip functional relationship 
is a dynamical equilibrium. 

At this point, we need a theory i) to support the idea 
that the freely evolving flow dynamics will effectively self- 
organize into a dynamical equilibrium state ii) to determine 
the q — ip relationship associated with this dynamical equi- 
librium iii) to select the dynamical equilibria that are likely 
to be observed. This is the goal and the achievement of 
equilibrium statistical mechanics theory, presented in the 
next subsection. 

b. The equilibrium statistical mechanics of Robert- Sommeria- 
Miller (RSM) 

The RSM statistical theory is introduced on a heuristic 
level in the following. There exists rigorous justifica tions 
of the theory, see e.g. iBouchet and Corvellec ( 2010f ) and 
references therein. 

A microscopic state is defined by its potential vortic- 
ity field g(r). If taken as an initial condition, such a fine 
grained field would evolve toward a state with filamenta- 
tion at smaller and smaller scales, while keeping in general 
a well defined large scale organization. Then, among all the 
possible fine grained states, an overwhelming number are 
characterized by these complicated small scale filamentary 
structures. This phenomenology gives a strong incentive 
for a mean-field approach, in which the flow is described 
at a coarse-grained level. 

For that purpose, and following iRobert and Sommeria 



(|199ll ). we introduce the probability p(a, r)da to measure 
a potential vorticity level a at a point r = (x,y). The 
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Table 1. Symbols and notations used in the text. 


OjIIlUOl 


n AnTllfl CY 11 
LLdllll L1UI1 


&x,y,z 


unit vectors in the meridional (%\ zonal (?/) and vertical (z) direction 


r = \ x i y) 


coordinate of a point, witli V7 == {dxi &y) 


i 
I 


L1111L. LUUI UlIldLL- 


it 


\-i c^o c? r\Y 7" nine i"Y i /H /^t/yy i i~yy o t i /Y y~i 
XYUooUj IdUlUo Ul UclUX HldLlUll 


fr, 
./() 


( ,nri a iq ty ci v fi in otov 

Vj'UI IUIIo LJdl dlllULd 


fl 
P 


r\ Qn pf" Q rir \7Y"\Y"T~1 1*17" 0" VCl H 1 PT1 T" 
yJl.CiLl.Vj Ldl y VUI LH_,ILV til dUICIIli 


D 


Hnmain win pro flip flow taVpt; nlHPP ( witli 7^ — 1 ^ 


n(r) 

U(l ) 






(fmp oTflinpn 1 "nnt"PTitifil \rctvi"if i 'i~i"\r npln 

\ llil^j pi. L (X11J.(_-VjI j UUlLjllLlCll VUI 11L11 V ±J.Lj±Vj1. 


(J 


1 Diro 1 r\i nArfin ct I i rm'f lpifir mn^n rv > — ^v~i 1 /~v~i 

1L.VC1 Ul JJOLCLlLildl VUI LjICILV, W1LI1 (J t ^-J — LXJ , ^L>J 


1-, 4> 


I c(~\ oyco n'r'imnrl 1 tnflcin "r\/~i"t"/~iTT^~i *i 1 \r/^T""t _ ir i iTT7' otiiH CTT'^dTnTiinptinn 
1 LUdlDt, dllltjU 1 IllUcxIl JJU UCLlLlcLl VUI lilLlUy cLIlU. o LI VjCHIIILUIC L1UI1 


Hi a 


Energy of the flow 


Y(er) 


global distribution of potential vorticity levels 




probability distribution function 




Mixing entropy 


j\ yj\ — iv 


IlUI llldjIlZiCll'lUIl CUllO LI diHIL 


u <y[P\ — l\ a ) 


constraint on global potential vorticity distribution 


fin] — W 


CUIlb LI ctlilL UII LIIL- Ulldty 




T , ci cfra n o"P mi i 1 ~i~ i nl i pt ticenn ci'f'pH ixri "t* n flip n r^ym tili r 7£ii~n~»Tn \l 
iJdgjI dli^L. UlLIlLipild dooU(_IclLt-U WILII Lilt- IIUI IIIdllZidLIUII ;V 


a(a) 


ijagiaii^e iiiuiLipnei associaieo wilii ine j^iouai voiLiciiy qisliiuulioii ±J a 


A 


\ nri* in o = ^a tyi i i t i i~\ 1 1 c\ y 1 cj c /"\ Piornr Tinr h t Iy n /in nvoi r 
Jjdgl dll^U IlllllUipild doDULldLtU WILII LIIL, v_,II(_,I g y O 




q - — ijj relation at equilibrium 




rcscalcd streamfunction 


TIM P 
J-[(p\, r 


Free energy functional and equilibrium free energy 


J VP) 


specinc 11 ee eneigy 


L/ — — Alt 


rcscaled Lagrange multiplier associated with, the energy constraint 


AT C A** A, 
1\1 — Jjj dvip 


constraint for (j) 




Lagrange parameter associated with M 


01: 02 


values of <p in a given phase 


?li 11 


Y 7"0 1 1 1 DO /"IT /Y 1 Yl O fYl Y r/YYl 1^\ t~\ O OZY 

values 01 c/ 111 d given pnaoe 




domain (and area) occupied by a given phase 


T 

lj 


TinnmAfpi' at t" Iy i"i lnfrirTOfri I ^, nfiiT'n n nnocAC 

pei inictei ui liil. 111 lci idee uclwccii piidoco 


T 


curvature radius of the interface 


V 


I q or ci n o'ri iyyi i 1 "t - i T - * 1 1 /yy" ci ccapi ci i~/Yri unth Trio c* f~\ ty cf - Y - ci i iy "t - / "i n /l „ urnriTi tyi i ti ityy i rzi yi <t / 
JjdKl all^L 1I1LU L1U11L.I dBOULldLLU WILII LIIU LUIlo LI dill L UII -ti^, ™ dlLll IllllllllllZiIIIg jLj 


T = Rt 


cnornin^tp fiprcm'i; tnr 1 lTitoT'iarr 1 

LUUI V.1J.1 LCAi V) Ci Vj L WOO Ui.lvj li.lLiVjl. L CL V_- v 


Y]et\< ) 


ipt Timfilp 

I U L/l Ulllb 


TJ. . — (A„ _ Ay-. )T? 
^ jet — Vv^2 s y iy JL 


velncitv of the iet 


Ti„ t = cBL(6o — S-\ "l 2 

J %nt ^i^i^\ t f f z 1/ 


Free energy of the interface, with c ~ 1 


C\ri\ — C l — ft 


pirnQ^vciiTil" r»n thp Iitiacit tyi cy tyi ptiti i tyi (? * in it"1 si 1 ■ r * titi a 1 1 

CUIIO LI dlllL UII LIIC IlllCdl IIIUIIICIIL LIIII It . LLIL LIdl , J Illldl ( 


tjf, V-pf 


latitude of the ring center of of jet latitude 


/? = /3/i? 2 


rescaled beta coefficient 




contribution of the beta term to the free energy 


Z(z) 


perturbation of the zonal interface 




First order corrections to the free energy 


in ij/ 


zonal and meridional extension of the closed domain 


-bring 


diameter of the ring 



probability density field p defines a macroscopic state of field, also referred to as coarse-grained, or mean- field, is 
the system. The corresponding averaged potential vorticity r 

5( r ) = / dcr <rp {<*,*) , (4) 
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with the average streamfunction ip defined by q = 'S/ 2 yj — 
ip/R 2 + (3y, and where £ =] — oo, +00 [. Many microscopic 
states q can be associated with a given macroscopic state 
rho. The cornerstone of the RSM statistical theory is the 
computation of the most probable state p eq , that maxi- 
mizes the mixing entropy given by the Boltzmann-Gibbs 
formula 

s Ip} = - dr / daplogp , (5) 
Jv Js 

while satisfying the constraints associated with each dy- 
namical invariant. The mixing entropy (O is a quantifica- 
tion of the number of microscopic states q corresponding 
to a given macroscopic state p. The state p eq is not only 
the most probable one: an overwhelming nu mber of mi- 
crost ates are effectively concentrated close to it (jMichel and 
1994). This gives the physical explanation and the predic- 
tion of the large scale organization of the flow. 

In the remaining of this paper, the term global entropy 
maximum or stable equilibrium state will be used for any 
global maximizer p of the entropy ([5]) satisfying the con- 
straints. The term local entropy maximum or metastable 
equilibrium state will be used for any state p that is a local 
maximizer of the entropy (O, satisfying the constraints. 

To compute statistical equilibria, the constraints must 
be expressed in term of the macroscopic state p: 

• The local normalization N [p] (r) = J E dcr p (cr, r) = 



equilibria characterized by a functional relation between 
potential vorticity and streamfunction: 



is d(T 



3 Aeri/>(r)-a((T) 



J s dcr e crA ^( r ) _Q ( ') 



(7) 



• The global potential vorticity distribution D a [p] = 
j v dvp{a,r) =7(0"), 

• The energy £ [p] = — \ J v dr J" s dcr p{a — (3y) tfj = E. 

Because of the overwhelming number of states with 
only small scale fluctuations around the mean field poten- 
tial vorticity, and because energy is a large scale quan- 
tity, contributions of these fluctuations to the total en- 
ergy are negligible with resp ect to the mean-field energy 
Robert and Sommerial ( 199l[ ). 

The first step toward computations of RSM equilibria is 
to find critical points p of the mixing entropy ([5]) . In order 
to take into account the constraints, one needs to introduce 
the Lagrange multipliers £(r), a(a), and A associated re- 
spectively with the local normalization, the conservation of 
the global vorticity distribution and of the energy. Critical 
points are solutions of: 



It can be shown that g a is a monotonic, increasing and 
bounded function of Xip for any global distribution 7(cr) 
and energy E. These critical points can either be entropy 
minima, saddle or maxima. To find statistical equilibria, 
one needs then to select the entropy maxima. 

At this point, two different approaches could be fol- 
lowed. The first one would be to consider a given small 
scale distribution 7(17) and energy E, and then to compute 
Rofitel statistical equilibria corresponding to these parame- 
ters. In practice, especially in the geophysical context, one 
does not have empirically access to the fine grained vor- 
ticity distribution, but rather to the q — ip relation ([7]) of 
the large scale flow. The second approach, followed in the 
remaining of this paper, is to study statistical equilibria 
corresponding to a given class of q — ip relations. 

More precisely, we will consider the class of q — if) rela- 
tions that admit an inflexion point, referred to as "tanh- 
like" relations, see figure [TJ When the global distribution 
7(cr) is a double delta function (a two level system), one 
can explicitly show that the q — ip relation is a tanh func- 
tion lBouchet and Sommeria ( 2002 ). but the actual class of 
initial conditions associated with ton/i-like relations is ex- 
pected to be much larger. A pragmatic point of view is that 
these tanh-\\ke relations are the one that allow for statisti- 
cal equilibria characterized by fronts of potential vorticity. 
In that respect, there are the relevant class of q— ip relations 
to describe either rings or zonal jets. 



V dp SS - X5S - / dcr a SD a - / dr (5N = , (6) 

JT, JV 

where first variations are taken with respect to p. This 
leads to p = Ncxp (Xatp (r) — a(cr)) where N is determined 
by the normalization constraint (J da p = l). Finally, us- 
ing Q, one finds that statistical equilibria are dynamical 



c. Simplification of the variational problem 

As explained in the introduction, a widely used method 
to solve a variational problems with many constraints (such 
as the RSM variational problem) is to consider a dual vari- 
ational problem, which has the same critical points as the 
initial one, but that is less constrained: any solution of the 
less constrained (easier to solve) probl em is a solution of 



the more constrained problem, see e.g. lEllis et al. I (l2000h . 



Let us for instance consider the ensemble of rescaled 
streamfunction fields <p = ip/R 2 , that satisfy the constraint 
J v dr <f> = M. This constraint is equivalent to the conser- 
vation of the first moment of potential vorticity J v dr q, at 
leading order in R. Let us then look for the minimizers of 
the free energy functional 



dr 



p 



^ (v^) 2 + / (0) - Hv 



(8) 



where /(</>) is a specific free energy to be defined precisely 
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in the next paragraph. This provides a variational problem 



F = min \ T [<t>] 

4> 



dr 6 = M 



(9) 



which is much simpler than the one of the RSM statistical 
theory, since only one constraint is kept. Critical points of 
this problem are solutions of ST — fi dr 8<j) = 0, for any 
perturbation 6(j>, where a is the Lagrange multiplier asso- 
ciated with the constraint. A part integration and the rela- 
tion q = R 2 V 2 cj>-<i)+(3y give ST = J dr (/'(</>) -<f>-q) 5<j>. 
Critical points satisfy therefore the relation q — f {<f>) — <f>— 
[i. These critical points are the same as the RSM critical 
points, given by equation ([7]), provided that 



f{cp) = g a {\R 2 4>) + d> + Li 



(10) 



This relation defines the specific free energy /. One can 
see that the ton/i-like shape of g a with a sufficiently steep 
slope at its inflection point leads to a double-well shape 
for /, as illustrated on figure [TJ This double- well shape is 
an important ingredient for the computation of statistical 
equilibria in the following. 



It has been proven by iBoucheti (|2008l ) that for each 
minimizer cj> of the variational problem (j9j), there exists a 
set of constraints E, 7(0-) such that <f> is the mean-field 
streamfunction of the RMS statistical equilibrium state p 
associated with the constraints E, j(a). In other words, 
any local (global) free energy minimizer cj> can be inter- 
preted as a local (global) entropy maximum state of the 
RMS theory. 

3. Potential vorticity homogenization and jets as 
statistical equilibria 

Assuming that there is no beta effect (/3 = 0), that f(4>) 
has a double- well sha pe, and considering th e limit ii<L 
(L is the domain size), Bouchet ( 2001 , [2008h found that the 
variational problem ((9]) becomes analogous to the Van der 
Waals - Cahn - Hilliard model that describes pha se sepa- 
ratio n and phase coexistence in thermodynamics (jModica 
1987). This formal analogy provides then an interesting 
physical interpretation of self-organization phenomena in 
geostrophic turbulence. 



a. Solutions of the Van der Waals 
problem 



Cahn - Hilliard variational 



At zeroth order in R, the function / (0) plays the dom- 
inant role in the free energy functional T given by ([H]). In 
order to minimize T, the streamfunction (f> must therefore 
be equal to one of the two minima of the specific free en- 
ergy /((/)) (the points 4>\ and </>2 on figured]). Each of these 
minima corresponds to one phase. Without the constraint 
Jp drc/> = M , one of the two uniform solutions <f> = <pi or 
4> = 4>2 would minimize T: the system would have only one 
phase. But in order to satisfy the constraint dr <j> = A/, 




Fig. 1. a) In this paper, we consider the class of q — <p re- 
lations (with 4> = 4>/ R 2 ) having a tanh-\ike shape, namely 
i)dccreasing with <f> ii) bounded for ip — ¥ ±00 iii) with a sin- 
gle inflection point. In addition, we assume that the slope 
at the inflection point is sufficiently steep, so that the q — <p 
relation represented as a plain line cross three times the 
dashed line q = —(f) — p. b) The double-well shape of 
the specific free energy / (4>) appearing in the expression 
of the free energy functional ((HJ. This function is related 
to the q — ip relation through equation (fT0|) . The Lagrange 
parameter (i is chosen such that the specific free energy 
of the two minima are the same (f(<pi) = f{4>2)), which 
allows phase coexistence. 

the system has to split into sub-domains: part of it with 
phase <j) — 4>i and part of it with phase <j) — 4>2- In terms 
of free energy minimization, the coexistence of these two 
phases is possible only if f(<f>i) = f{4>2)- Using equation 
(fT0|) . one can always choose the Lagrange parameter p to 
satisfy this condition, sec figure [1] for a graphical interpre- 
tation. In physical space, the area occupied by each of 
the phases is denoted Ai and A2 respectively, see figure 
[2] These values are fixed by the constraint J v dr <j> = AI, 
which gives at leading order (j)iAi + (^2^2 = M and by the 
geometrical constraint Ai + A2 = 1 (where 1 is the area of 
the domain). 

The interface between the sub-domains characterized 
by 4>i and <p2 corresponds to an abrupt variation of stream- 
function. The term R 2 (V</>) 2 in the expression ([8]) of the 
free energy J- is negligible everywhere except around this 
interface, on a typical width of order R. The interface 



G 



a) Initial condition b) Formation of two phases c) Interface minimization 




Fig. 2. Resolution of the Van der Waals variational problem ©. a) Example of an initial condition for the potential 
vorticity field. Note that this initial condition could as well have many levels of potential vorticity. b) At zeroth order 
in R, <f> takes two values 0i, <p2 on two sub-domains A\ and A2 corresponding to the coexistence of two phases of 
homogenized potential vorticity qx, q2- These sub-domains are separated by strong jets of typical width lj et = R and 
velocity Uj e t = (<f>2 — <t>i)R- c) The actual shape of the structure, or equivalently the position of the jets, is obtained by 
minimizing its perimeter L for a fixed value of A2. 



is therefore associated with a strong and localized jet di- 
rected along this interface, with a typical velocity Ujet = 
(4>i — fa) R, and a typical width R. The actual jet profile 
is computed in Appendix A, by minimizing the free energy 
associated with this profile. The jet gives always a positive 
contribution to the free energy: 



H) RL , withc~l 



(11) 



In order to minimize this intcrfacial free energy, the 
perimeter of the jet L must be minimal, taking into account 
the constraints given by the fixed areas A\ and A2 . We thus 
look for the curve with the minimal length that bounds a 
given surface. The solution of this classical problem is that 
the interface is either a circle or a straight line. 

To conclude, the computation of the statistical equi- 
libria predicts i) the formation of two phases of constant 
streamfunction, with strong and localized jets at these in- 
terface ii) the velocity profile across of these jets iii) the 
shape of the interface, which is determined by an isoperi- 
metrical problem: the minimization of the interface length 
for a fixed enclosed area. 

b. Link with potential vorticity homogenization theories 

The tanh-likc shape of the q — ip relations imply that 
sub-domains of constant streamfunction arc also sub-domains 
of constant coarse grained potential vorticity. It means 
that the potential vorticity is homogenized in each sub- 
domain. Statistical mechanics provides therefore a phys- 
ic al explanation for the pot ential homogenization theory 
of Rhines and Younel (1982), without invoking any dissi- 
pation mechanism. 

In the case of freely evolving 1.5 layer QG dynamics, 
statistical mechanics predicts not only the spontaneous for- 
mation of regions where potential vorticity is homogenized 



at a coarse grained level, but also the shape of the inter- 
face between these regions, corresponding to jets, where 
vorticity gradients arc confined. 

Statistical mechanics arguments also account for the 
irreversible nature of mixing: an overwhelming number of 
fine grained microscale are associated with a given coarse 
grained equilibrium state. The only effect of a weak small 
scale dissipation process would be to smooth out locally 
fine-grained fluctuations of potential vorticity, leaving un- 
changed its coarse-grained structure. 

Note that the formation of two subdomains of homog- 
enized coarse-grained potential vorticity is essential to en- 
sure the energy conservation. Importantly, any eddy pa- 
rameterization based on local down-gradient diffusivities 
could not represent this homogenization process: it would 
lead to the formation of a single phase of homogenized 
potential vorticity, which would in general not satisfy the 
energy constraint. 

These results are complementary to previous work fo- 
cusing on the dynamics of potential vorticity mi xing, using 
chaotic advection theory (jPierrehumberti 1 1 9 9 lh . Chaotic 
advection theory has been proven successful to account for 
many observed feature of potential vorticity mixing, but 
unlike statistical mechanics, it provides in general no pre- 
diction for the final state of the large scale flow, especially 
when there is no scale separation between mean and eddies 
in the initial condition. 

4. Application to oceanic rings 

Observations show that mesoscale oceanic rings exist 
everywhere in the ocean, and particularly near western 
boundary currents, where high levels of eddy kinetic en- 
ergy levels are la rgely associated w ith their prese n ce, se e 
e.g. lOlsonl (Il99lh for a review and iMorrow et al.l (J2004J) ; 
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Chelton et al. (l2007h for recent altimetry measurements. 
Both cyclonic and anticyclonic rings propagate westward 
at speed (3R 2 . They also present an additional small merid- 
ional drift, poleward for cyclonic rings and equatorward for 
anti-cyclonic rings. 

Am ong others , co ntributions of McWilliams and Flierll 



1979).lNoJ (|l981l) . of lFlierll (|1987I) and of lCushman-Roisin et 



1990|) have provided insights on the dynamics and on the 



mechanisms responsible for the generation of the rings, in 
a large class of models and configurations. Various dynam- 
ical mechanisms accounting for the formation of the rings 
have been pointed out: vortex shedding above topography, 
pinch off process during the non-linear evolution of a mean- 
dering jet, boundary layer separation, or large scale orga- 
nization of initially small turbulent disturbances. Despite 
these very different generation mechanisms, the striking re- 
semblance between observed oceanic rings in very different 
regions of the oc ean, which h as long been recognized as a 
surprising result (|01sonlll991l ). suggests that at least some 
aspect of these coherent structures can be studied indepen- 
dently of their generation mechanism. This is precisely the 
interest of statistical mechanics, which accounts for spon- 
taneous formation of circular structures surrounded by a 
jet, i.e. the self-organization of mesoscale turbulence into 
rings. 

a. The westward drift of the rings 

To show that westward propagating circular rings can 
be interpreted as equilibrium states, two important ingre- 
dients must be taken into account: i) the beta effect (3y 
ii) a domain invariant by translation in the zonal direc- 
tion. A drifting ring could not exist as statistical equilibria 
in a closed domain, since it would be destroyed when ar- 
riving on the western boundary. The zonal translational 
invariancc of the problem has important consequences. It 
is shown in Appendix B that a change of Galilean reference 
frame in the zonal direction translates as a beta effect in 
the expression of potential vorticity. Moreover, in a refer- 
ence frame moving at velocity —j3R 2 e Xl the beta effect is 
exactly canceled out. We conclude that in a domain invari- 
ant by translation in the zonal direction, statistical equi- 
libria obtained by the minimization of the Van-Dcr-Waals 
Cahn Hilliard variational problem © without beta effect 
are also statistical equilibria with beta effect, but drifting 
westward at speed V ~ —(3R 2 . 

b. The poleward drift of cyclones and the equatorward drift of 
anticyclones 

If the flow actually reaches a local statistical equilib- 
rium, then not only the ring is composed of an homoge- 
nized region of potential vorticity, but also the background 
flow. In figure [3] is represented the case of an isolated patch 
of potential vorticity (q = qi within the ring of area Ai cen- 



tered on y = 0) on a beta plane (q = (3y elsewhere). This 
situation is common in the ocean, for instance when Agul- 
has rings arrive in quiescent regions of the Atlantic ocean. 
Since the background potential vorticity is not homoge- 
nized, this state is not a statistical equilibrium state. 

It is shown in the following that the observed asymmet- 
ala ic small meridional drift of cyclonic and anticyclonic rings 
can be understood as a tendenc y for the system to reach the 



can be understood as a tendency tor trie system to reach the 
statistical equilibrium (see also lSchecter and Dubinl ( 2001 ') 
for a similar argument in the context of beta plane turbu- 
lence). One needs for that purpose to consider the conscr 
vation of the linear momentum: 

C = / dr qy , 



which is the dynamical invariant associated with the zonal 
translational symmetry. The linear momentum of the ini- 
tial condition is, at lowest order in R: 

£<«/?/" dr y 2 - (3 f dr y 2 , 
Jv Ja z 

where Ai is the initial area of the ring. Considering the 
limit of small rings compared to the domain size, the first 
term of the right hand side dominates the second one, and 
C l ps (3L x Ly/2A. Assuming the statistical equilibrium is 
reached in the final state, the flow is made of two phases 
of homogenized potential vorticity: the background phase, 
with value q = qi,, and the ring's phase of area Af, centered 
at latitude y = Vf, with value q = qf (same sign as q{) and 
with | qf | > |<jt>|. The linear momentum of this final state 
is 



C f ps q b / dr y - (q f - q b ) / dr 
Jv J A f 



y 



Finally, the linear momentum conservation C l = & gives 



24 v 



(if - lb) y/Af . 



Physically, this means that for statistical mechanics reason, 
the background potential vorticity has to be homogenized, 
which leads to a loss of linear momentum that must be com- 
pensated by a latitude shift of the ring center. Rings with 
negative potential vorticity ((qf — qt) < 0) decrease their 
latitude (yf < 0), while rings with positive potential vor- 
ticity (qf — qi,) > increase their latitude (yf > 0). Then, 
in order to reach the statistical equilibrium, the ring has 
to drift northward if it is made of an initial positive poten- 
tial vorticity patch, and southward if it is made of an initial 
negative potential vorticity patch. This corresponds always 
to a poleward drift for cyclonic structures and an equator- 
ward drift for for anticyclonic structur es, just as what is re- 
ported from altimet ry measurements (jMorrow et al.l 12004 : 
Chelton et al~ll2007l) . 



8 




Fig. 3. Explanation of the meridional drift of rings, as a tendency to reach the statistical equilibria. On the left, initial 
conditions, with a white disk of positive potential vorticity in the upper panel, and a black disk of negative potential 
vorticity in the lower panel. In both cases the disk evolves on an initial beta plane with no background flow. On the 
right, the corresponding statistical equilibrium. 



c. Conclusion: oceanic rings are local statistical equilibria 

In the ocean, the scale separation between the size of 
the rings and the Rossby radius of deformation is satis- 
fied only to a limited extent. This scale separation has 
been assumed for technical reasons only: it allows for ex- 
plicit analytical computations of the equilibria. The results 
obtained in this limit actually apply for far more general 
situations. This is confirmed by numerical computations of 
RSM equilibria, in which rings are obtained as local equi- 
libria even when the scale separation is not satisfied, see 
e.g. figure |U 

To conclude, the quasi-circular shape of oceanic rings 
and their westward propagation suggest that these coher- 
ent structures can be interpreted as local statistical equilib- 
ria. The existence of a meridional drift shows a departure 
from the prediction of the equilibrium theory. However, 
the fact that this drift can be interpreted as a tendency 



to reach to equilibrium state shows that these structures 
remain close to an equilibrium state. 

Rings are local (metastable), and not global statisti- 
cal equilibria of the equivalent barotropic model: a global 
equilibria would imply the coalescence of all existing rings 
into a single large scale vortex, in order to minimize the to- 
tal interface between the different regions of homogenized 
potential vorticity. 

5. Application to mid-basin eastward jets 

Another region of the ocean where strong jets of typ- 
ical width given by the Rossby radius of deformation R 
are localized along an interface separating two regions of 
homogenized potential vorticity is the inertial part of mid- 
latitude eastward jets, as the Kuroshio or the Gulf stream 
currents. The inertial part of these currents is located 
in the regions where the western boundary currents sep- 
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Stream function Velocity modulus 




Fig. 4. Circular vortex as a statistical equilibrium of the 
QG model, with R < L r i ng . Although analytical compu- 
tations are carried in the limit R -C Lringi the results are 
expected to hold when this scale separation does not exit. 
It is a circular patch of (homogenized) potential vorticity in 
a background of homogenized potential vorticity, with two 
different values. The velocity field (right panel) has a ring 
structure. The width of the jet surrounding the ring has 
the order of magnitude of the Rossby radius of deformation 
R. 

arate from the coastline and self-organize downstream into 
a strong eastward jet. Because mid-basin eastward jets 
fill a large part of oceanic basins, and because the exis- 
tence of a western boundary is an essential ingredient for 
their formation, one must look for statistical equilibria in 
a close domain in that case. In view of the applications 
to mid-basin ocean jets, we assume a situation in which 
the global distribution of potential vorticity is symmetric: 
two phases characterized by symmetric values of stream- 
function (4>i = —<p2 = Ujet/ (2R)), each of them filling half 
the domain area {A\ = A2 = 1/2). We ask in this sec- 
tion whether configurations with mid-basin eastward jets 
arc statistical equilibria. 

a. Without beta effect, mid basin eastward jets are statistical 
equilibria of the QG model 

The value — (pi t 2 for the two coexisting phases is not 
compatible with the boundary condition = 0. As a con- 
sequence, there exists a boundary jet in order to match a 
uniform phase = 0x,2 to the boundary conditions. Just 
like interior jets, treated in section [3J these jets contribute 
to the first order free energy, which gives the boundary 
jet structure and shape. The symmetry of the problem 
(0i = — 02 ) implies that boundary jets of each phase give 
the same contribution to the free energy. Because the 
boundary length is a fixed quantity, the free energy min- 
imization amounts to the minimization of the interior jet 
length only, just as in previous subsections. The interior 
jet position and shape is thus given by the minimization of 
the interior jet length with fixed area A2 = 1/2. 




Fig. 6. Strcamfunction of the solution presenting an east- 
ward jet with beta effect (red: positive values, blue: neg- 
ative values) , associated with case (a) of figure [5j The jet 
width is of order R. This solution is a statistical equilib- 
rium for L x < iriU-jet/P) 1 / 2 . 

The jet has to be straight or circular. There are three 
possible interface configurations with straight or circular 
jets: i) the zonal jet configuration (jet along the x axis) 
with L — L x , ii) the meridional jet configuration (jet along 
the y axis) with L = L y , iii) and an interior circular vortex, 
with L = \Z2tt . The case i) of a zonal jet is a global 
interface minimum (and then a global equilibrium state) 
if and only if the aspect ratio satisfies L x /L y < 1. For 
L x /L y > 1, these solutions become metastable states (local 
entropy maximum). 

We conclude that without beta effect, mid-latitude east- 
ward jets are statistical equilibria. Because of the symme- 
try 0i = —02 , solutions presenting eastward and westward 
jets are equivalent: westward jets are also statistical equi- 
libria. 

b. With beta effect, eastward jets becomes metastable or unsta- 
ble 

Contrary to the case of the zonal channel, the beta 
effect can not be cancelled out by a change of Galilean 
reference frame in the case of a closed domain. One can 
therefore not avoid taking into account this term in the 
computation of the equilibrium free energy. One can read- 
ily see on the expression (jSJ of the free energy, that when 
j3 0, the additional term Fp = —j3 j v <ir <f>y breaks the 
symmetry ±q. The westward jet case (with < on the 
southern part of the domain and > on the northern 
part) is more favorable in terms of free energy minimiza- 
tion than the eastward jet case (with > in the southern 
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Fig. 5. a) Eastward jet configuration b) Westward jet configuration c) Perturbation of the interface for the eastward 
jet configuration, to determine when this solution is a local equilibrium. Without beta effect, both the eastward and the 
westward configurations are entropy maxima. With positive beta effect the westward jet becomes the the global entropy 
maximum, and the eastward jet becomes metastable provided that beta is small enough. 



part of the domain and <p < in the northern part): west- 
ward jets become the only global equilibria for (3 > and 
aspect ratio L x /L y > 1. 

Let us be more precise by considering the limit of small 
beta effect, with the scaling (3 ~ R(3. With that scaling, 
the equivalent topography does not play any role at zc- 
roth order in the variational problem ((9]). We thus still 
conclude that phase separation occurs, with sub-domains 
of fixed areas A\ and A 2 , separated by jets whose trans- 
verse structure is described in Appendix A. It is shown 
in this same Appendix that the interface gives a contribu- 
tion Tint = cRL (fa — fa) 2 ■ Using the zeroth order result 
4> = fa on sub-domain A\ and <f> = fa on sub-domain 
A21 one obtains also J-p = — (fa — fa) J A dr y , plus an 
unimportant constant. Finally, the total first order contri- 
bution to the free energy is 



Tr = cR I 



fa) 13 / dr y 

J An 



(12) 



Recalling that first variations of the length are propor- 
tiona l to the inverse of the curv ature radius r of the inter- 
face ( Gelfand and Fomhl 1967 ) , the minimization of (fT2j) , 
with fixed area A2 gives 



(fa-fa)Rpy + i] 



cR{ 



l) 2 



(13) 



where 77 is a Lagrange parameter associated with the con- 
servation of the area A 2 . 

We conclude that zonal jets (i.e. an interface at latitude 
V = Vjet, with infinite curvature radius r) are solutions to 
this equation for rj = —R (fa — fa) f3yj et . This shows that 
eastward and westward jets described in the previous sec- 
tion are therefore still critical points of entropy maximiza- 
tion. 



The eastward jet configuration is the one with the re- 
gion A 2 below the line y = at the center of the do- 
main. To determine if this configuration is a local sta- 
tistical equilibria, let us consider perturbations of this in- 
terface (given by the line y = l(x)), while keeping constant 
the area occupied by both phases, see figure [S] There are 
two contributions competing with each other in the ex- 
pression (|12p of the free energy Tr. Any perturbation in- 
creases the jet length L = J dx yl + (I') 2 , where I' = 
dl/dx, and then increases the first term of the free energy 
02) by ST mt = cR(fa - fa) 2 J dx (l'f. Any perturba- 
tion decreases the second term of the free energy (fT2)) by 
STp = —R (fa — fa) (3 J dx I 2 . If the eastward jet solution 
is not a free energy minimum, it exists a perturbation of 
the interface leading to negative variations of the free en- 
ergy STr = STint + 5Fp. Let us consider the particular 
case I = lkSui(kTTx/L x ) where k > 1 is an integer. Using 
Uj e t — (fa — fa) R and (3 = Rf3, the condition 8Fr < 
gives f3 > cl/jet (kn/L x ) . The most unfavorable case is for 
the smallest value of k 2 , i.e. k 2 = 1. 

It leads to the necessary condition (3 > cUj et TT 2 / L 2 for 
the eastward jet solution to be unstable (in term of sta- 
tistical mechanics). It can actually be shown, using less 
straightforward considerations, that it is also a sufficient 
condition for instability. The destabilizing effect of in- 
creasing values of (3 contrasts with its stabilizing e ffect in 
classic al criteria for barotropic instability, see e.g. IVallid 
( 2006h . It has actually been shown that such eastward jet 



solutions can simultaneously be unstable for statis tical me- 
chanics, and stable for n on-linear perturbations (|Venaillel 
20081 : IVenaille et al.ll201lh . 



For a fixed value of /3, eastward jets are local free en- 
ergy maxima if the domain zonal extension is smaller than 

1/2 

a critical value: L x < n (cllj e t/ f3) . The strcamfunction 
of such a state is presented in figure [6] For jets like the 
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Gulf Stream, Ujet ~ 1 w.s -1 and /3 
Using c ~ 1, the critical domain length scale upon which 
eastward jet become unstable is L x « 300 km. This length 
is smaller than the typical zonal extension of the inertial 
part of the Kuroshio or Gulf Stream currents, but not by an 
order of magnitude, which suggests that these structures 
are marginally unstable. The instability is consistent with 
the fact that strong meanders and pinch off process occurs 
downstream of oceanic eastward jets. But the marginal 
nature of this instability is also consistent with the overall 
robustness of the global structure of the flow, which be- 
comes a statistical equilibrium when the extension of the 
jet is small enough. 

6. Conclusion and prospects 

The aim of this paper was to present a point of view 
complementary to existing approaches that deal with co- 
herent structures in the ocean. It was shown that the 
RSM statistical mechanics provides a unified framework 
that may be useful to study mesoscalc and basin scale in- 
ertial flows. 

The theory gives a physical explanation and a predic- 
tion for the self-organization of large scale oceanic coher- 
ent structure, independently of the underlying generation 
mechanism. It predicts the formation of subdomains of 
homogenized potential vorticity, with intense jets at the 
interface. Mesoscale rings can be interpreted as local equi- 
librium states of the RSM theory. Their shape and their 
drift can be understood in this framework. Mid-basin east- 
ward jets are found marginally unstable states of the RSM 
theory, consistently with observations of these jets. 

The interest of this approach relies on its generality (it 
docs not depends on a particular flow configuration) and 
on its ability to describe qualitatively different observed 
regimes of self-organization, such as rings and zonal jets. 
The present study was achieved in the framework of a 1.5 
layer QG model, which is too simplistic to describe oceanic 
eddies quantitatively; however, generalizations and further 
investigations in the framework of more complex models 
can be built upon these results. 

A caveat of this approach is that forcing and dissipation 
are not taken into account in the framework of the equi- 
librium theory: the input of the RSM theory is given by 
the dynamical invariants. In the case of mesoscale rings, 
even if the dynamics can be considered close to an equi- 
librium state, forcing and dissipation play an important 
role in setting these dynamical invariants. In the case of 
basin scale jets, their marginal instability suggest that one 
can not avoid taking into account forcing and dissipation 
mechanisms to explain these structures. So far the inertial 
part of wind driven circulation has been mostly studied 
from the point of view of bifurcation theory, starting from 
a highly dissipated ocean and decreasing progressively fric- 



tional parameters, see e.g. iDiikstra and Ghill (|2005l ). We 
argue that this problem can be tackled with another point 
of view, starting from the purely inertial limit (this pa- 
per), and a dding small forcing and dissipa tion (future work 
built upon Bouchet and Simonnet ( 20091 )). These two ap- 
proaches are complementary and may be combined in the 
future in a more comprehensive non-equilibrium theory. 
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APPENDIX A 



Computation of the jet profile 

At leading order, minimization of the free energy leads 
to the formation of subdomains of constant streamfunction 
<f> = 4>x and (f> = <p2- The interface between these subdo- 
mains are associated with strong and localized jets. Let us 
assume that the curvature radius of the interface is much 
larger than i?, which allows us to neglect what happens 
along the interface at leading order. Calling r — Rr the 
coordinate in a direction along the normal to the interface, 
the jet profile <pj et (r) across the interface must be such 
that it minimizes its contribution to the total free energy 
JSJ). The jet profile is therefore determined by solving a 
one dimensional variational problem: 



Tim = LR min 







d(j)jet \ 




1 J — oc 




dr J 



f(4>jet) 



(Al) 

where L is the perimeter of the jet and T n t the free 
energy associated with the existence of this interfacial jet. 
Critical points of this variational problem are states that 
cancel the first variations of the free energy with respect 
to (pint- They are solutions of 



d 2 <p 



dr 2 



df 



(A2) 



Making an analogy with mechanics, if <pj et would be 
a particle position, r would be the time, equation (| A2|) 
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would describe the conservative motion of the particle in a 
potential — /. In order to connect the two different phases 
in the bulk, on each side of the interface, one has to consider 
solutions with boundary conditions (f> — > <f>i for r — > — oo 
and </>—>■ 4>2 for r — > +00. It exists a unique trajectory 
with such limit conditions. In the particle analogy, it is 
the trajectory connecting the two unstable fixed points <\>\ 
and </>2, corresponding to the two bumps of the potential 
— / (see figured]). 

The energy (dcpjet/dr) 2 /2 — f{4>jet) is conserved dur- 
ing the evolution of <p with time r. Using this conser- 
vation property and the boundary condition <j) — > 02 for 
r — > +00, one obtains f(4>j e t) = {d4>jet / dr) 2 /2, plus an 
unimportant constant. Injecting this expression into the 
variational problem (|A1[) . one obtains 

/+00 
(d^et/dr) 2 dr = cRL (0 2 - (j)i) 2 , 
-00 

with c ~ 1. An important physical consequence is that the 
jet at the interface always gives a positive contribution to 
the free energy of the equilibrium state, which is of order 
R and proportional to the interface length L. 

APPENDIX B 

Galilean invariance and beta effect 

In the case of a zonal channel, the QG equations (fTJ) 
are invariant over a Galilean transformation in the zonal 
direction 

x — > x' = x — Vt , y — > y = y , t — > t = t . 

The velocity is transformed as v — > v' = v — Ve x , which, 
using the relation v = e z x V^, gives the transformation 
law for the streamfunction %j) — > tj)' = tj) — Vy. From the 
expression q = Atp — ip/R 2 + [3y, one obtains finally the 
transformation law for the potential vorticity q — > q' = 
q + Vy/R 2 . Thus the expression for the dynamics in the 
new reference frame is 

da' 

+ v' • VY = 0, with v' = e 2 x V V , 

of 

and g ' = V 'V-|i+(^+^)y'- 
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